Performance of Linear Field Reconstruction 
Techniques with Noise and Uncertain Sensor 

Locations 

Alessandro Nordio*, Carla-Fabiana Chiasserini, Emanuele Viterbo 
Dipartimento di Elettronica, Politecnico di Torino 
C. Duca degli Abmzzi 24, 1-10129 Torino, Italy 
Phone: +39 0115644183, Fax: +39 0115644099 E-mail: 

{alessandro .nordio, carla. chiasserini, emanuele . viterbo} Spoilt o . it 

* Corresponding author 



o 
o 



in 

q 

O 



> 

p 

o 
o 



X 



Abstract — We consider a wireless sensor network, sampling a 
bandlimited fleld, described by a limited number of harmonics. 
Sensor nodes are irregularly deployed over the area of interest 
or subject to random motion; in addition sensors measurements 
are affected by noise. Our goal is to obtain a high quality 
reconstruction of the field, with the mean square error (MSE) of 
the estimate as performance metric. In particular, we analytically 
derive the performance of several reconstruction/estimation tech- 
niques based on linear filtering. For each technique, we obtain 
the MSE, as well as its asymptotic expression in the case where 
the field number of harmonics and the number of sensors grow 
to infinity, while their ratio is kept constant. Through numerical 
simulations, we show the validity of the asymptotic analysis, even 
for a small number of sensors. We provide some novel guidelines 
for the design of sensor networks when many parameters, such 
as fleld bandwidth, number of sensors, reconstruction quality, 
sensor motion characteristics, and noise level of the measures, 
have to be traded off. 

EDICS: DSP-RECO Signal reconstruction, DSP-SAMP 
Sampling, SEN-APPL Applications of sensor networks, SEN- 
FUSE Data fusion from multiple sensors, SPC-PERF Perfor- 
mance analysis and bounds 



1. Introduction 

Wireless sensor networks are often used for applications 
like environmental and traffic control, habitat monitoring, or 
weather forecasts [1], which require to sample a physical 
phenomenon over an area of interest (the sensor field). In this 
paper, we consider a set of sensors communicating with a sink 
node, through either single- or multi-hop communications. 
Each sensor locally samples the physical field, while the sink 
collecting all samples is in charge of reconstructing the signal 
of interest. 

We assume that initially sensors are either located at pre- 
defined positions, or, if randomly deployed over the network 
area, their location can be estimated at the sink node (see [2]- 
[4] for a description of node location methods in sensor 
networks). We do not deal with spatio-temporal correlation, 
but consider a fixed time instant and focus on the spatial sam- 
pling and reconstruction of the sensor field. We note that, in 
general, sensors provide an irregular sampling of the observed 
phenomenon. This may be due to various reasons: random 
deployment of the nodes, environment characteristics that 
bias the network deployment, sensors entering a sleep mode, 
inaccuracy in sensor positioning, or nodes movement [5]. In 
all these cases the sink has to reconstruct the field from a 
collection of samples that are irregularly spaced, different 
from the classical equally (or regularly) spaced sampling. 

The problem of signal reconstruction from irregular sam- 
ples has been widely addressed in signal processing, where 
several efficient and fast algorithms have been proposed to 
numerically reconstruct or approximate a signal [6], [7]. The 
problem we address in this work, however, is different; the 
questions we pose are: 

(i) How do noisy measures and inaccurate knowl- 
edge of the sensor positions affect the quality of the 
reconstructed signal? 

( ii) How can we trade off system parameters like 
measurement noise, field bandwidth, signal recon- 
struction quality and number of sensors? 

To answer these questions we analyze two different models 
of the monitoring system that account for the quality of the 
measurements performed by the sensors and differ in the 



accuracy with which the sensor positions are known at the 
sink node. More specifically, the model denoted as Model 
A refers to the case where sensors are fixed, the sink has 
perfect knowledge of the sensor positions, but the sensor 
measurements are affected by error In the second model, 
named Model B, besides noisy measurements, we consider 
that the sensors position varies around an average value, and 
only the average location of the nodes is known at the sink. 
Examples where this model applies are observation systems 
using surface buoys [8], underwater robots located at different 
depths [9], [10], dropsondes or low-cost unmanned platforms, 
as in [11]. 

For each of these models, we use as field reconstruction 
techniques some linear filters that are commonly employed 
in signal detection and estimation, and we evaluate the mean 
square error of the resulting estimate. 

We find that a key parameter for the network performance 
is the ratio /3 of the field number of harmonics to the number 
of sampling sensors. In particular, there exists a value of 
this ratio, beyond which the performance of all considered 
reconstruction strategies degrade significantly, even for low 
values of noise level and limited uncertainty on the sensor 
positions. To obtain an acceptable reconstruction quality when 
j3 is large (i.e., the number of available sensors is limited 
compared to the field bandwidth), reconstruction techniques 
that exploit some knowledge of the measurement noise and of 
the jitter in the sensors position must be employed. 

We also carry out an asymptotic analysis of the system as 
the field number of harmonics and the number of sensors grow 
to infinity, while their ratio (3 is kept constant, and we show 
that this is an effective tool to study the system performance, 
even when the number of sensors is small. Finally, we find 
a lower bound to the mean square error that can be achieved 
by any of the considered techniques, both under Model A and 
Model B. 

The remainder of the paper is organized as follows. In 
Section |ll] we present our assumptions and the system models 
under study. Section |lll] highlights our contribution with re- 
spect to previous work. Section |IV] introduces the performance 
metrics and provides some mathematical tools necessary for 
our study. Model A and B are analyzed in Sections |V] and [Vl] 
respectively. Finally, in Sections IVIII and IVIIII we summarize 
our main results and draw some conclusions. 

II. Assumptions and System Models 

Let us consider a one-dimensional bandlimited field s{x) 
represented by 2M + 1 harmonics as 



s{x) 



^/2M 



1 ^ 



(1) 



k=-M 



The field is observed within one period interval [0,1) and 
sampled by r sensors placed at position^ x = [xi, . . . , Xr]^ , 

'Column vectors are denoted by bold lowercase letters, matrices are denoted 
by bold upper case letters. The (k, q) entry of the matrix X is denoted by 
(X)ii;q. The n X n identity matrix is denoted by In, the generic identity 
matrix is denoted by I, and the conjugate transpose operator is denoted by 
(■)t 



Xq E [0,1), q = l,...,r which are in general not equally 
spaced. The signal samples are denoted by the column vector 
s = [s{xi), . . . , s{xr)f^- The field discrete spectrum is given 
by the 2M + 1 complex vector a = [a_A/, • • • , ao, • • • , om]"^- 
The complex numbers represent amplitudes and phases of 
the harmonics in s{x). We can think of M as the approximate 
one-sided bandwidth of the field. 

We assume that the entries of x are i.i.d. uniformly dis- 
tributed random variables in [0, 1). The extension to a multi- 
dimensional field can be easily obtained, as discussed later in 
this section. 

We define (3 as the ratio of the number of harmonics 
which describe the field to the number of sensors, i.e., /3 = 
(2Af + l)/r. This is an important parameter in our analysis. 
Note that the number of sensors r also corresponds to the 
sampling rate; thus, the number /? is the ratio of twice the field 
bandwidth to the sampling rate (frequency). In particular, in 
regular sampling theory, exact reconstruction is achieved for 
/? € [0,1) and, if a Nyquist regular sampling interval were 
used, we would have: /3 = 1. 

We consider AI to be known, and the random vector a to 
have zero mean and covariance matrix E[aa^] = o'^I2a/+i, 
where cr^ corresponds to the field average power spectral 
density. 

The value s of the field at positions x depends on the 
spectrum a through the expression 



Gla 



(2) 



where Gx is the {2M + 1) x r generalized Fourier matrix 
defined as: 



(Gx)fcg 
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k = -Af, 
g= 1,... 



(3) 



V2M+T 

The dependence of the matrix Gx on the position vector x 
is clearly indicated by its subscript. When the samples are 
equally spaced in the interval [0, 1), the matrix ^//3Gx is a 
unitary matrix (i.e., PGy^G]^ ~ l2M+i)- The above system 
model refers to a uni-dimensional field where sensor positions 
are determined by a scalar variable. However, the extension to 
the multi-dimensional case can still be easily obtained since 
the relation between field spectrum and samples in a band- 
limited multi-dimensional field can be expressed in a matrix 
form similar to Q, where only the structure of the matrix Gx 
differs. 

Finally, we assume that sensor field measures are sent to 
a processing unit, the so-called siyik node, whose task is to 
provide an estimate of the sensed field. Since we focus on the 
reconstruction of the physical field, we consider that sensor 
transmissions always reach successfully the sink nod^ 

By relying on the assumptions discussed above, we study 
the following two systems. 

« Model A: Fixed sensors, perfect knowledge of the 
sensor positions, noisy measures 
In this model, sensors have a fixed position, given by 
the vector x and known at the sink node, but each 

-Note that this is a fair assumption since, when ARQ or EEC techniques 
are used, the information either is correctly retrieved at the sink or it is lost. 
The latter case coiTesponds to reduced value of r 



sensor provides a measure of the field affected by additive 
noise with zero mean and variance cr^ [12]. The additive 
noise approximates the errors affecting the measurement 
procedure [13]. 

The measures vector can therefore be written as: 

p = s + n = G],a + n (4) 

where s is the true field and the zero mean noise vector 
is denoted by n, with covariance matrix E[nn^^] = (^fj-r- 
• Model B: Sensors with jittered positions and noisy 
measures 

In this case each sensor moves around an average position 
Xq (q = 1, . . . ,r), i.e., the sensor positions are given by: 
X = X + ^, where: E[x] = x and S is the displacement 
of the sensors with respect to their average location 
X. Note that our problem differs from the well known 
problem of jittered sampling (see e.g., [28]), since we 
deal with irregular sample locations. The displacements 
6q, q = l,...,r, are modeled as independent zero 
mean Gaussian random variables with variance ag and 
E[<^5^] = (Tg Ir- For convenience and neglecting the edge 
effects, we consider mod {xq, 1) so that Xq falls in the 
observation interval [0, 1). The vector p of measures is 
still given again by (|4|. Also, noise, displacement, and 
field spectrum are assumed to be uncorrected, hence 
E[n^t] = E[nat] = E[a^^] = 0, and the sink has perfect 
knowledge of x. 

III. Our contribution with respect to previous 

WORK 

Given a network where sensors can enter a low-power 
operational state (i.e., a sleep mode), the work in [14] presents 
an algorithm to determine which sensor subsets should be se- 
lected to acquire data from an area of interest and which nodes 
should remain inactive to save energy. A similar problem is 
addressed in [15], where an adaptive sampling is described, 
which allows the central data-collector to vary the number of 
active sensors, i.e., samples, according to the desired resolution 
level. The optimal sensor density that minimizes the network 
energy consumption, subject to constraints on the quality of the 
reconstructed signal and network lifetime, is studied in [16]. 
Note that in our work we consider an irregular topology, which 
may be caused by nodes moving into a sleep state; however 
we do not directly address energy efficiency or scheduling of 
the node sleep/activity periods. 

In [17], the authors consider a uni-dimensional field, uni- 
formly sampled at the Nyquist frequency by low-precision 
sensors. The impact on the field reconstruction accuracy of 
quantization errors and node density is evaluated. The effect of 
random error sources affecting the ADC, besides quantization, 
is investigated in [13]. In our work we consider an additive 
noise that models errors due to the measurement procedure 
as well as errors due to the ADC, but we do not specifically 
focus on the latter issue. 

The impact of medium access control (MAC) protocols on 
the reconstruction of a signal field is investigated in [18]. Both 
deterministic and random MAC schemes are considered, and 



performance are derived as the number of received packets 
and the experienced SNR vary. 

Related to our work is also the literature on spectral analy- 
sis [30], [31], which deals with the problem of recovering the 
amplitude of sine waves immersed in noise. Note, however, 
that techniques such as MUSIC do not estimate phases; thus, 
we do not compare with such techniques since our linear 
filtering reconstruction yields the estimate of both amplitudes 
and phases. 

The field reconstruction at the sink node with spatial and 
temporal correlation among sensor measures is studied in 
[12], [19]-[21]. In particular, in [21] the observed field is 
a discrete vector of target positions and sensor observations 
are dependent. By modeling the sensor network as a channel 
encoder and exploiting some concepts from coding theory, 
the network capacity, defined as the maximum ratio of target 
positions to number of sensors, is studied as a function of 
noise, sensing function and sensor connections. The paper by 
Dong and Tong [22] focuses on signal reconstruction from 
possibly random samples, as we do. However, two major 
issues make our work significantly different from [22]. Dong 
and Tong indeed assume that the exact sensors locations 
are known and that the central controller always receives a 
sufficiently large number of samples. These assumptions allow 
an interpolation method, which is used in Dong and Tong's 
work, to provide good performances. In our case, instead, 
even in the asymptotic analysis, the ratio of the number of 
harmonics to the number of samples is kept constant and, 
hence, interpolation may be highly inefficient as we will show 
in the following. 

The problem of reconstructing a band-limited signal from 
an irregular set of samples at unknown locations is addressed 
in [23]. There the signal is oversampled by irregularly spaced 
sensors; sensor positions are unknown but always equal to an 
integer multiple of the sampling interval. Different solution 
methods are proposed, and the conditions for which there 
exist multiple solutions or a unique solution are discussed. 
Differently from [23], we assume that the sink can either 
acquire or estimate the sensor locations and that sensors are 
randomly deployed over a finite interval. 

Finally, in our previous work [24] some conditions on 
the irregular topology of the sensor network are identified, 
which allow for a successful signal reconstruction, both under 
deterministic and random node deployment. In particular, in 
[24] the spectrum estimate a, computed by the sink, is obtained 
by applying to s the Moore-Penrose pseudo-inverse of the 
matrix Gx, i.e., a ~ (G^G^) GxS. The system model 
adopted in [24] is ideal in the sense that the reconstruction 
algorithm has perfect knowledge of the vector x and neglects 
noisy measures: the failure in reconstruction (i.e. a 7^ a) is 
only due to the bad conditioning of the matrix GxG^^. in 
relation to the finite machine precision. In this work, instead, 
we propose to apply Unear filters to the field reconstruction 
and consider the following causes of quality degradation: (i) 
noisy measures, and ( ii) uncertainty at the sink on the sensors 
position. 



IV. Preliminaries 

Here we describe the techniques we use for field recon- 
struction, and define the performance metrics employed for 
assessing the effectiveness of these techniques on the quality of 
the reconstructed field. Finally, we provide some mathematical 
tools necessary for the analysis of the models under study. 

A. Reconstruction techniques 

Several reconstruction techniques have been proposed in the 
literature, which amount to the solution of a linear system 
(see [6], [7] and the references cited therein). A widely used 
technique consists in processing the measures p by means of 
a Unear filter, B, which is an r x {2M + 1) matrix and is a 
function of the system parameters known at the sink. In this 
case, the estimate of the field spectrum is given by: 



Btt 



(5) 



The system model in (|4]i is similar the one employed in 
multiuser communications [25] or multiple antennas commu- 
nications [26], [27]. In those cases p is the received signal, 
the matrix Gx plays the role of spreading matrix or channel 
matrix, a is the transmitted signal and n is the channel noise. 
By relying on the results obtained in those fields, for each 
system model we propose and compare some reconstruction 
techniques characterized by different matrices B: the matched 
filter (MF), the zero forcing (ZF) filter and some linear filters 
minimizing the mean square error (LMMSE) [25]. In the 
field of multiuser detection, the MF simply correlates the 
received signal with the desired user's time reversed spreading 
waveform, thus it does not take into account any other users 
in the system or channel dynamics. The ZF filter counteracts 
multiuser interference but it ignores the presence of channel 
noise. The LMMSE solution minimizes the squared error 
between the received and transmitted signals, thus accounting 
also for the channel noise; it becomes the zero-forcing solution 
when no noise is present. Note that the advantage of the MF 
with respect to the ZF and LMMSE filters is that no matrix 
inversion is needed; while, between the ZF and the LMMSE 
filter, clearly the best performance in terms of minimum square 
error is given by the LMMSE, however the advantage of the 
ZF filter is that it does not require any knowledge of the noise 
component (see [25] for further details). 

B. Performance metrics 

Given the spectrum estimate Q, the field can be recon- 



structed as: 
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j27Tkx 



k=-M 



As a measure of the quality of the reconstruction, we consider 
the mean square error (MSE) of the estimate of s{x), which 
is given by: 



MSE = E 



\s{x) — s(.t)P d.T 



We observe that computing MSE as above is equivalent to 
computing E[||a — a|p]. Indeed, we have: 



\s{x) — s{x)\^ Ax 
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Therefore, in the following, for a given vector of sensor 
positions x, we consider the MSE defined as: 



MSEx 



E |_!|a-a||- 
2M + 1 



2i\/ + 1 



^Tr{*x} 



where 



*x - — E[(a-a)(a-a)t] 



(6) 



(7) 



is a (2M+1) X (2M+1) matrix, the operator E[-] averages with 
respect to all random variables of the model, and Tr{ } is the 
trace operator. Also, in Q we exploited the fact that, for any 
vector V, we have: E [||v|p] = E [Tr{vvt}] Tr {e [vv^] }. 

Next, we consider the vector x to be random. In this case 
a more appropriate performance metric is the average MSE, 
normalized to cr^, i.e.. 



EfMSE, 



MSEa 



where MSEx is as in ^ and E[ ] averages over the realizations 

X 

of X. 

When the parameters M and r grow to infinity while 
the ratio /3 = {2M + l)/r is kept constant, we define the 
asymptotic average MSE as: 



MSEo 



lim MSEa 

A/,r — >-\-oo 

2M + 1 



(8) 



Our results will show later that MSEoo gives a very good 
approximation of MSEat) already for small values of M. This 
is a common feature of asymptotic analysis based on random 
matrices [26]. 



C. Some mathematical tools 

1 ) The fimctional (f): Let us first consider an nx n Hermitian 
random matrix X and the functional: 

0(X)= lim -E[Tr{X}] 

Using O and (|8]l, the asymptotic MSE can be written as: 

MSEoo = 0(*x) (9) 



In our analysis we use the following results on the functional 
4>{-). First, we notice that: </)(!) = 1. Secondly, we can prove 
that, if g{x) is an analytic function defined in x > ther|l 



0(.9(X))=E[.9(e)] 



(10) 



where ^ is a random variable with the asymptotic eigenvalue 
distribution of X. The proof is given in Appendix U 

2) A simple expression for G^+s- As will be clear in 
Section |VI] in the analysis of Model B many parameters are 
functions of the matrix Gx, where x = x + (5. It is thus useful 
to derive an expression of Gx as a function of Gx, in order 
to separate the random part (5 of x from the constant part x. 
From (O, the (fc, q) entry of Gx is defined as: 



(Gx)fcg - 
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V2AI + 1 



^— j27rA,'2;g 
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A useful expression of Gx in terms of Gx is given below. 

Lemma IV.l For any vector x of size r, let the {k,q) entry 
of the matrix Gx be 



V2M+T 



for k ~ — Af, . . . , M, and q ~ 1, . . . , r. Let the size r column 
vectors x, x, and S be such that x = x + 5, then 



oo ^ 

Gx = V ^W"GxA" 



(11) 



n=0 



where A = diag(<5) is an r x r diagonal matrix, and W is a 
(2M + 1) X (2M + 1) diagonal matrix with {'W)kk = —'ftuk. 

Proof: The proof is given in Appendix |ll] 

V. Analysis of Model A 

Here we consider the case where sensor positions are fixed 
and known at the sink but the field estimates are degraded by 
noisy measures. We analyze three different linear filters: the 
matched filter, the zero forcing filter and the minimum mean 
square error filter [25]. In all cases, for any fixed x, the filter 
matrix B is deterministic. Thus, using (|5]l, and Q we 
obtain: 



- E [||Bt(Gta + n)-a||2] 
(B^Gi - I)(GxB - I) + aB^B 



(12) 



where 



SNR,„ = - = ^ 



is the signal-to-noise ratio on the measure. The MSB expres- 
sion specialized to the different filters is given below. 



A. Matched filter 

As a first solution, we choose B as the filter matched to Gx. 
The MF is optimal when the collected samples are equally 
spaced, that is when the rows of y^SGx are orthonormal 
vectors and v^Gx is a unitary matrix (i.e., /3GxG^ = 
l2j\/+i)- Thus, we choose: 



Bt = /3Gx 



(13) 



Recall that Gx depends on the position vector x that, under 
Model A, coincides with the actual sensor positions. Indeed, 
in the absence of noise and for equally spaced sensors, we 
have the spectrum estimates perfectly match a, i.e., 

a = Btp = /3GtGxa = a 

By replacing ( fT3] l in ( fT2] i. we obtain the following expres- 
sion for 'S'x: 



*x = Rx 



I 



(a/3 - 2)/3Rx (14) 

where Rx = GxG^.. From the definition in (|9]l, the asymptotic 
MSB, averaged over the random vector x, is given by: 

MSEoo - </'(*x) 

- /320(R2) + 0(I) + (a/3-2)/30(Rx) 

Notice that the second term on the right hand side reduces to 
1 since </)(I) = 1. Applying (fTOl i. first with g{x) = x^ and 
then with g{x) = x, we obtain: 



MSEoo = E[A2] + 1 + (a/3 - 2) E[A] 



(15) 



where A > is the random variable with probability density 
function (pdf) f\^p{x), distributed as the asymptotic eigenval- 
ues of /3Rx. In [24] it is shown that, for any positive integer 
p, E[Ap] is a polynomial in /3 of degree p — 1. In particular 
E[A] = 1 and ^[\^] = 1 + /3. We therefore obtain: 

MSEoo = /3(a + 1) (16) 

B. ZF filter 

The expression of the ZF filter for the system in dUi is: 

Bt ^ R-iGx (17) 

Notice that, by its definition, the ZF filter does not exploit any 
information on the noise contribution (such as a"^). However, 
this reconstruction technique takes into account the fact that 
the collected samples are not equally spaced and, hence, that 
-y//3Gx is not a unitary matrix. 

By using ( fTTI i in ( fT2b . the matrix '^y^ becomes: 



(18) 



Using the definition in Q and applying ( fTot with g{x) = x^^, 
the asymptotic MSB, averaged over the random vector x, can 
be written as: 



MSEoo = a(/)(Rx^) = a/3 E 



(19) 



^ Note the small abuse of notation when using g(-) for both scalar and 
matrix argument 



We can make the following observations on the behavior of 
the MSEoo: 



1) Since 1/A is a convex function, then E[l/A] > 1/E[A]. 
In [24] it is shown that E[A]=1, thus it resuhs: MSEqo > 
a(3. 

2) We have: MSEoo a(3 < +oo only for (3 e 
[0, P"), with w 0.35. Indeed 



E 



"1' 




A 


-I 



-fx,p{x) dx 



In [24] it has been empirically observed through Monte- 
Carlo simulation that for x <C 1: 

where the exponent a(/3) is a decreasing function of f3 
for /3 e [0, 1], and a(/3) = 1 for P = 13". Given that, for 
any positive constant c, we have: 



where the integral in the right hand side (and there- 
fore ([T9T l) does not diverge if and only if a{(3) > 1, that 
is P < (3* . This observation gives us a fundamental limit 
to the minimum number of sensors required to perform 
reliable reconstruction with the ZF filter. 



C. LMMSE linear filter 

A more efficient solution is to employ the filter B that 
provides the minimum MSB (LMMSE). By assuming that the 
signal-to-noise ratio SNR„i is known to the sink and exploiting 
this information for the filter design, the expression of the 
LMMSE filter [25] for Model A in © is given by: 



Bt = (R, + aI)-iGx 



(20) 



We highlight that this reconstruction technique accounts for 
both the fact that the collected samples are non-uniformly 
spaced and the presence of the measurement noise. 
Substituting (l20t in (fT2l i. we obtain: 

*x = a (Rx + (21) 
Using ([Tot with g{x) = {x + the asymptotic MSE is: 



MSE. 



E 



a(3 



A + a(3 



Note that: 



E 



a/3 



Also note that E[a/3/(A 



> 



a[3 



E[A + a/3] 1 + al3 



(22) 



(23) 



a/3)] < 1, since A > 0. Given that 
the LMMSE filter provides the minimum MSE, from (l23T l it 
turns out that, for a given /3, a/3/(l + a/3) is a lower bound 
for the performance of all linear reconstruction techniques. 
We summarize the main results of this section in Table U 



TABLE I 

Results obtained under Model A 
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D. Results 

In Figure [T] we compare the average MSE obtained using 
the MF, ZF, and LMMSE filters, when (3 varies and a = 1/2 
(i.e., SNR,„ = 3dB). The points labeled by "MSEav MF", 
"MSEav ZF" and "MSEav LMMSE" have been obtained 
generating 100 realizations of the measures (|4]i with M ~ 40, 
computing the estimates as in (|5]) and averaging the square 
error ||a — a|p. These points are superimposed to the solid 
curves labeled by "MSEoo", representing the asymptotic MSE 
and obtained evaluating ( fTst . ( fT9t and (l22l) . respectively. 
Notice that computing closed form expressions for E[l/A] 
in ( fT9] l and E[a/3/(A + a/3)] in ( |22] | is still an open problem 
since a closed form expression of the distribution of A is 
unknown. Thus, for a given /3, the value of these asymptotic 
expressions have been obtained pseudo-analytically, averaging 
over the eigenvalues A obtained by several realizations of 
the matrix /3Rx, with il/ = 200 which yields a very good 
approximation of the asymptotic case (see [24]). 

We observe an excellent agreement between the asymptotic 
analysis and the numerical results; this shows the validity of 
the asymptotic analysis even for values of M as low as M = 
40. We also note that, for both the filters, higher values of 
MSE are obtained as (3 increases. Finally, the LMMSE filter 
provides the best performance, while the MSE of the ZF filter 
shows a vertical asymptote for /3 = /3*, in agreement with the 
closed form analysifl 
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Fig. 3. MSB obtained through the LMMSB filter, plotted versus SNR„i, for 
/3 = 0.1, 0.2, 0.4, 0.6, 0.8 and M = 10 



Figure |2] shows the MSB versus SNR,„, for (i = 0.2. The 
behavior of the asymptotic MSE is represented by the curves 
labeled by "MSEoo" while the average MSE obtained through 
numerical analysis is denoted by the label "MSEav "■ The 
curves have been obtained using the same procedure as for the 
results in Figure [H using M = 10 for MSEau and M = 10 
for MSEoo computation. Again, note the tight match between 
analytical and numerical results. For all techniques, the MSE 
decreases as the SNR„i increases. The MF however provides 
very poor performance, even for high SNR^. In particular, as 
SNR,„ tends to infinity, it shows a horizontal asymptote with 
MSEoo = /?. 

Besides linear filtering, another technique for estimating 
the spectrum a is based on interpolation [22]. The idea is to 
interpolate the measures p to a regular sampling grid defined 
by the vector x' where = (g — l)/r, q = 1, . . . , r. The 
interpolated vector p' is then multiplied by the matrix /3Gx' ■ 
Notice that in this case v^Gx' is unitary i.e. /^Gx'G]^., = I, 
since x' represents an equally spaced sampling. In the figure 
the dashed line labeled "Linear interp." shows the performance 



obtained using linear interpolation. The MSE has a horizontal 
asymptote for high SNR„i. While it outperforms the MF, it 
clearly shows poor performance for high SNR, compared to 
ZF and LMMSE techniques. 

Figure |3] presents the performance of the LMMSE filter 
obtained evaluating (|22] | for different values of /?, as the SNR,„ 
varies. In agreement with the results presented in Figure [T] the 
performance of the LMMSE filter degrades as /? increases, 
while, as expected, it improves as the SNR™ increases. 

Example 1: We need to estimate the number of sensors 
required to sample a field with M — 100 harmonics. Each 
sensor provides samples with SNR^ = 30 dB. 
We choose to employ the LMMSE filter, which provides 
the best performance. Looking at Figure [3] if we allow 
an MSEav of 3 • 10"^, then we need (3 = 0.4, i.e., r = 
(2M + l)//3 sa 500 sensors. By doubling the number of 
sensors (/3 = 0.2), MSEav drops to 3 • 10"**. 



VI. Analysis ob Modbl B 

Here we consider the case of sensors with jittered positions 
and average position, x, known at the sink node. The true 
sensor location is: x = x + ^, where (5 is a random vector, as 
defined in Section HI] The reconstruction algorithm employs 
the matrix B, which is a function of the known average 
positions x. For any given x and B, similarly to the MSE 
becomes: 



MSEi 



E [||. 

a.n.o 



2Af + 1 



2Af + 1 



Tr{*i} (24) 



where 



4 E [(a-a)(a-a)t] 
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Bt ( E[GiG, 



al B 



ifBtB] 



25R <! E[G, 
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(25) 



where 3?{ } represents the real part of the argument. 

To proceed further we need to compute the averages over 

the displacements S, i.e., we need the expression of E[Gx] 

s 

and E [G^^-Gx] as functions of Gx, whose derivation is given 

s _!_ 
in Appendix Hill We have: 



E[Gx] = CGx 
s 



and 



E [Gt Gx] = Gt C^Gi 
s 



1 - 



Tr{C^} 

2M + 1 



I 



(26) 



(27) 



where C is a {2M + 1) x {2M + 1) diagonal matrix with 
(C)fcfc = C5(-j27rfc), k = -M,...,A/, where Cs{-) is 
the characteristic function of the displacements. Under the 
assumption that 6 has a zero mean Gaussian distribution we 
have (C)fcfc = cxp{-27r'^Pa^s), k = -M, ...,M. 



Using ( |26] | and dZTI i in (IZST i we obtain: 

*i = Bt (Gtc2Gi + 7l) B-23fi{CGxB} + I 

(28) 

where7 = l + a + ^igl. 

In the following, in the case of the LMMSE filteo we 
first consider that the variance <jg of the sensor movement is 
unknown at the sink and, hence, the sink assumes the sensors 
to be fixed (i.e., S = 0), while running the reconstruction 
algorithm. Then, we consider that is known and the 
reconstruction algorithm employs a filter that exploits such 
an information to minimize the MSE (this case is referred to 
as "LMMSE for known af\) 

Finally, we remark that, while in Model A the filters used 
for signal reconstruction are functions of the matrix Gx with 
X known to the sink, in Model B only the mean value of 
the sensor positions x is known and, hence, the filters are 
computed using Gx instead of Gx. 

A. Matched filter 

If the sink node employs the MP in ( fTSl ) as function of x 
(i.e., Bt = /3Gx), then, using ( l28t , we obtain: 

*i = /g^R^C^Ri + 7/32Ri _ {CRx} + I (29) 

This result holds for strictly positive a^. Note that, for aj — O 
(no sensor motion), we have C = I and j = a; thus ( l29l l 
reduces to (fT4l i. 

Equation ( l29b refers to the MSE obtained with a given 
vector x; we are now interested in deriving the asymptotic 
expression for the MSE. Note that ( |29] l is a function of both 
Rx and C, and contains terms of the form C^gCR.^) with 
g{x) = l,x,x^ and p = 0, 1,2; also the matrix Rx depends 
on M and r, while the matrix C depends on M and aj. The 
definition of the asymptotic MSE in ^ refers to the case 
where the number of harmonics AI and the number of sensors 
r grow to infinity with constant ratio (3; if this is directly 
applied to ( |29] l. information losses may arise. Indeed, we have: 
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(30) 



and thus all terms depending on the matrix C would vanish 
regardless of the value of a^. On the contrary, in a realistic 
situation we expect to obtain high reconstruction quality when 
the standard deviation of the motion (as) is smaller than 
or comparable to the average sensor separation (1/r), and a 
significant degradation of the reconstruction quality when as is 
much larger than the average sensor separation. To distinguish 
such different conditions, we define the signal-to-noise ratio 
on the motion as: 

{l/rf 1 



^Recall that the MF and ZF techniques, by their definition, do not require 
any information on cr^ and (j| 



where lo = asr. We then redefine the asymptotic MSE as 
the limit of the average MSE for A/, r +oo, with constant 
/3 = {2M + l)/r and constant oj = asr. In this case. 
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where iy{x) = /Acrfl^TTx) / (nx). Notice that i/{0) 



(31) 
1 and 



</)(7) = l + , 



0. Also, we have: 
Tr{C2} 



lim 

M,r^+oo 2M 
I3.UJ 



1 



1 + a - v{l3uj) (32) 



Using the new definition and 
of the MSE becomes: 

MSE^ = Tr{*x} 

= 0(/32R-C2Ri- 



, the asymptotic expression 



7/32Ri-2/33fi{CRi} + l) 



= (C2r|)+/320(7R.)_ 2/30 (CRx) + 1 
= /32</,(c2)</,(R?) + /32</>(^)0(R^) 

-2/3(/.(C)0(Rx) + 1 
= z^(/3w) E[A'] + /90(7) E[A] 
{(iuj/^) E[A] + 1 

-2iy{f3u/V2) + 1 
= P{l + a) + i'{l3uj)-2iy{Puj/V2) + l (33) 

Here we used the following facts: 

• (/)(3fi{CRx}) = 0(CRx) since Rx is Hermitian and C 
is real and diagonal; 

« Tr{XiX2} = Tr{X2Xi} for any square matrix Xi and 
X2; 

• (/)(CPR?) — (?!)(CP)0(R|) for any positive integer p 
and q. This assumption holds only if C and Rx are 
asymptotically free [26]. Since asymptotical freeness is 
in general very hard to prove, we will simply verify the 
validity of such assumption through numerical results. 

. E[A2] = 1 + /? and E[A] = 1 (see [24]); 
Equation ( |33] ) reduces to ( fTSI l for a; = 0, while it reduces to 
MSEoo = 1 + /3(1 + a) for uj = +00. 

B. ZF filter 

In this case the sink node employs the ZF filter in (1% but, 
knowing only the average value of the sensor positions, the 
filter results to be a function of x: B^ = Rr^Gx, and the 
matrix 'S'x can be written as: 



*x = ' + (c - 1)' 



(34) 



We observe that, when ag ~ (no sensor motion), we have 
C = I and 7 = a, thus ( ijib reduces to ( fTsT l. 



Using OTT i and ( |32] |. the asymptotic MSE is; 
MSEoo = 0(7Ri' + (C-I)2) 

A 

+ 1 + I'iPuj) - 2iy{Puj/V2) (35) 

Equation (l35T l reduces to ( fT9] l for a; = 0, while it reduces to 

MSEoo = 1 + /3(1 + Q') E[l/A] for u ^ +oo. 

C. LMMSE filter neglecting ag 

If the sink employs the filter in ( |20] | computed using x (i.e., 
= Ar^Gx, where = + aT), then the matrix *x 
in (l28T i becomes: 

*x = A7iRi(c2Rx + 7l)A7i „ 25R{CRxA^i} +1 (36) 



For cr| = (i.e., C = I and 7 = a). 



reduces to dJTJ. 



Using the properties described in Section IVI-AI the asymp- 
totic MSE is: 

' A2 - 



MSE. 



1 



(A + a/3)2 
A 



(A + ap)''_ 

Equation (IJTT i reduces to (l22l i for = 0, while it becomes: 

MSEoo = 1 + + a) E[A/(A + a^f] for +cx). 

D. LMMSE filter for known cr| 

We now consider the linear LMMSE filter optimized for 
the case where (t| is known at the sink. We find the optimal 
B minimizing Tr{\I/i}; that is, we null the derivative of (|24] | 
with respect to B. We employ the following properties that 
hold for any square matrix X [29]: 



A5RTr{XB} = Xt 
ATr{BtXB}=2XB 



if X = xi' 







Then, we have: 

9MSEx(B) _ (gIc^Gx + 7I) B Sagct C 
dB 2M + 1 2M + 1 

Solving for B, we obtain the expression of the LMMSE filter 

Bt = (CRiC + 7!)"^ CGi (38) 

Substituting ( l38T l into (l28l l. we have: 

*x = 7(CRiC + 7l)-i (39) 

In this case an explicit expression of MSEqo is hard to obtain. 
However, we were able to find the following lower bound that 
turns out to be very tight, as shown by the results presented 
in the following section 



MSE. 
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+ 0(7) 
a — v{(3uj) 
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(40) 



where to derive the last expression we exploited OTI) . (l32b . 
( [Tol l and the fact that E[A] = 1. 
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Fig. 4. Performance of the ZF filter for /3 = 0.2 and M = 10, when cr^ is 
neglected 




Fig. 5. Performance of the LMMSE filter (20) versus SNR,„, for P = 0.2 
and M = 10, when (t| is neglected 



E. Results 

We now show the performance of the filters analyzed under 
Model B. Regarding the ZF filter ( fTTI i. Figure |4] compares 
the asymptotic MSE evaluated through ( |35] | (represented by 
solid lines and labeled by "MSEoo") against the average MSE 
(represented by points and labeled by "MSEav")- The MSEav 
is obtained by generating 100 realizations of the measures as 
in (|4]i, with M = 10, computing the estimates as in (O and 
averaging the square error ||a — ajj^. The MSE is shown in 
the log scale plotted versus SNR„i, for j3 = 0.2 and different 
values of SNR^^;. 

Similarly, Figure |5]presents the performance of the LMMSE 
filter ( l20t . Here the curves labeled by "MSEqo", generated 
through evaluation of ( |37] |. and the points in the plot, labeled 
by "MSEav", have been obtained as for Figure |4] 
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Fig. 7. Performance comparison of the LMMSE filter neglecting o-f (|20) 
against the LMMSE filter with perfect knowledge of (t| , as SNRm varies and 
for = 0.2 and M = 10 

In both the plots the solid line labeled by "SNR^r = +00" 
refers to the case where = 0, i.e. ^ = 0, and correspond 
to the performance provided by Model A under the same 
conditions. 

The excellent match between the asymptotic results and the 
numerical simulation confirms the validity of the asymptotic 
analysis as an effective tool to characterize the performance 
of the reconstruction techniques. 

Also, comparing Figures |4] and |5] we observe that the 
performances of the ZF and the LMMSE filters are similar 
for SNR,„ > 10 dB for any value of SNR^., while, for lower 
SNR„,, the LMMSE filter outperforms the ZF filter. 

Figure |6] compares the performance of the LMMSE fil- 
ter ( |38] |. which has knowledge of aj, with its lower bound ( l40l i 
(dashed lines), as SNR,„ varies. We consider (3 = 0.2 and 
different values of SNR^^. Notice that the lower bound is very 
tight, especially for high values of SNR^^. The points in the 
plot, labeled by "MSEav" have been obtained as for Figure |4] 
using M = 10. Here, as well as in Figure |5] the the case 



SNRj: = 00 (solid line) is shown, and corresponds to the 
performance of the LMMSE filter for signal model A. Indeed, 
for SNR2: = 00 (i.e., cr| = and 5 = 0), we have C = I and 
7 = a, and ( |39] | simplifies to (|2T]) . 

Figure [7] compares the performance of the LMMSE filter 
labeled by "LMMSE" (solid lines), and of the LMMSE 
filter (I38]l, labeled by "LMMSE erf' (dashed Hnes), for the 
same parameter setting as in Figure |6l For the considered 
value of (3 (/3=0.2), the filter in dSSl ) outperforms the sim- 
pler filter ( l20t for any value of SNR„i and SNR^,, but the 
performance gain is always negligible. 

Figure [8] shows the performance of the LMMSE filter ( |20] | 
neglecting erf obtained through evaluation of dsTI i for SNRj^ = 
lOdB (dashed Hnes) and SNR^, = 20 dB (soHd Hnes), and for 
(3 = 0.1, 0.2, 0.4. While the MSE^o of the LMMSE filter ^ 
always tends to 1 for small values of SNR^ (i.e., large values 
of a), for high SNR^ (i.e., low a) its behavior depends on /?. 
Indeed the term E[A/(A-f q;/?)^)] on the right hand side of dSTT i 
reduces to E[l/A] for a ^ 0. As explained in Section IV-BI 
E[l/A] diverges for /3 > /3* w 0.35 and so the MSE (see the 
lines with o markers in the plot). This behavior is more evident 
as (3 increases and the MSE is large, for any SNRm- These 
results, however, are of no interest from the application point 
of view since a system characterized by such poor performance 
is not working. 

Finally, Figure |9] compares the performance of the LMMSE 
filter ^ labeled by "LMMSE erf' (dashed lines) and the 
performance of the LMMSE filter (EQIi, labeled by "LMMSE" 
(solid lines), for SNR^ = lOdB and /3 = 0.1,0.2,0.4. 

In general the filter dSST l performs always better than fil- 
ter (I20I 1. In particular, for (3 < (3* the two filters show very 
similar performance, while, when (3 > (3*, the filter dSST l does 
not diverge for high SNR,„. This is shown in Figure |9] where, 
for P > 0.35 and high values of SNR^, the advantage of 
exploiting the knowledge of (t| becomes evident. 
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Fig. 9. Performance compaiison of the LMMSE filter neglecting (jg against 
the LMMSE filter with perfect knowledge of (t|, as SNRm varies, for 
SNR^ = 10 dB, and M = 10 



Example 2: Consider r ~ 1000 buoys deployed in water 
and equipped with sensors, which provide noisy measures 
with SNRm = SOfiffi. Buoys are moving but the variance 
(j| = 10"'' is unknown to the reconstruction algorithm. We 
need to estimate the maximum number of harmonics of the 
field that the network can sample and reconstruct with an 
average MSE lower than 5 10^'^. 

Since SNR„j is known to the reconstruction algorithm 
while ag is not, we employ the LMMSE filter in Sec. lVl-Cl 
We have: SNR^ = l/(cr|r2) = 10. Looking at Figure [H 
we notice that, for SNR,„ = 30 dB, values of MSEav 
lower than 5 10"'^ can be obtained only for (3 < 0.2. The 
maximum number of harmonics is then M = (?'/3— 1)/2 w 
100. 



Example 3: Consider a network of sensor with jittered 
positions characterized by (3=0.4 and SNK^ = 10 dB, and 
assume that these values are known to the reconstruction 
algorithm. We want to determine which type of sensor 
devices should be used in order to minimize the MSEav 
In other words, we ask ourselves how accurate the sensor 
measurements need to be (clearly, more expensive devices 
provide a higher SNR,„j. 

Since SNRa;=10dB is known to the reconstruction algo- 
rithm, we can employ the LMMSE filter given in ( |38] |. 
Looking at Figure |9] we notice that the performance of 
the filter for (3 — 0.4 shows a horizontal asymptote 
corresponding to an average MSE of 5 • 10~^. Thus, an 
SNR„i = 25 dB is enough to achieve the best performance. 



VII. Summary of Results 

Our main results for the system models A and B are as 
follows. 

Model A (fixed sensors and noisy measures): 



• for a given /?, the MSE provided by any of 
the reconstruction techniques is lower bounded 
by a/3/(l + a(3) and worsen with increasing (3 
(i.e., the ratio of the number of harmonics to the 
number of sampling sensors); the MF in ( fTST l 
is the only filter which does not require matrix 
inversion, however it provides poor performance 
in all of the considered cases; 

• the ZF filter provides high quality performance 
only for high SNR,„ (namely, SNR„ > 10 dB) 
and (3 < 0.35; 

• the performance of the LMMSE filter, instead, 
is good moderate values of SNR„i and f3 < 1. 

Model B (sensors with jittered positions and noisy mea- 
sures): 

• for a given (3 the MSE provided by any of the 
reconstruction techniques is lower bounded by 

• the performance of all reconstruction techniques 
worsen with increasing (3 and SNR^; 

• the advantage of exploiting the knowledge of 
SNR:c in the filter design is negligible for low f3 
and low SNR™, while it is of fundamental im- 
portance to obtain a high quality reconstruction 
for P > 0.35 and large values of SNR™. 

VIII. Conclusions 

We addressed the problem of reconstructing band-limited 
fields from measurements taken by irregularly deployed sen- 
sors, and we studied the effects of noisy measures and jittered 
sensors positions on the reconstruction quality. We analytically 
derived the performance of several linear filters in terms of the 
MSE of the field estimates. We also studied the asymptotic 
MSE, obtained as the number of harmonics and the number 
of sensors grow to infinity while their ratio (3 is kept constant. 
We found that the asymptotic analysis is an effective tool to 
characterize the performance of the reconstruction techniques 
even for a small number of sensors, and we investigated the 
impact that the parameter f3 has on the system performance. 
In [24] we observed that random sampling without any type 
of noise would require more than twice the sampling rate 
{(3 < 0.5) of minimum regular sampling {f3 = 1) to get 
a reliable reconstruction (without ill conditioning problems) 
with high probability. The number of sensors further increases 
{(3 < 0.2) compared to regular sampling when measurement 
noise (model A) and sensors position jitter (model B) are 
present. 
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Appendix I 
Proof of (fTob 



Let X = AA^ be an Hermitian matrix where A is an 
arbitrary random matrix. Let g(-) be an analytic function, 
defined for x G M.~^ that can be written as: g{x) ~ J2t=^ cLix"^, 
with finite coefficients a^. Considering that (/'(X'') is the p-th 
moment of the asymptotic eigenvalue distribution of X, i.e., 
^(X*') = E['C'^] where ^ is the random variable distributed 
as the asymptotic eigenvalues of X, and the continuity of the 
function we have: 



0(5(X)) 



^a,^(X') = ^a,E[f] 



+ 00 



+ CXD 



E 



+ 00 



Li=0 



E[5(0] 



Appendix II 
Proof of Lemma ITvTI 

Using (fTTT i the fcg-th entry of Gx is 

oo ^ 

= E;^EE(w")..(G.).,(A) 

n— h j 

oo ^ 

= E;;i(w"W(G^).,(A") 



I L. 

ri=0 

exp((W)fcfc(A),,)(Gx)fe9 

exp {-'fiMq) {Git)kq 

^^^^ _^ ^ CXp {~j2TTkSq) exp {-]2TTkXq) 

1 



exp (— j27rfca;q) 



x/2l7TT 
which matches its definition. 



Appendix III 

"""^'t^Bjvfi^U TATION OF E [Gx] AND E [G^-Gx] AS FUNCTIONS 

OF Gi 

We derive here the expressions of E [Gx] and E [G]j.Gx] 
as functions of G<r. 



a) Computation of E[G:x]- Using Lemma IIV.T] we have and 



E [Gx] = 



OO _j oo ^ 

ni ' Til 



.ji=0 



Tl = 



.An) (m) 

[GtG.] = ^^^^Z 



7 'rrj ' 



The average of A" is given by: E [A"] 



n=0 
m=0 



r X r 



x" fs{x) dxir ~ /i'-"^ Ir-, where I,, is the 
identity matrix and z^*^") is the n-th moment of S. Hence, 

EfGJ = y_W"Gi/ x"fs{x)dx 

V^W"/5(a;)dxGx 
■■^ — ' n! 

exp(xW)/5(a;)dxGi = CGi 

-OO 



E /, ^ diag(Z) (41) 



n=0 
m=0 



The first term of ( |4TI ) yields 
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The (fc, fc)-th element of Z is given by 



cxp{x'W)fs{x) dx 



(z; 



kk — 



is a (2M + 1) X (2M + 1) diagonal matrix and 
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kk 



cxp(— j27rA:a;)/i(x) dx 



E(Gt)fe„(Wt«W™)„,(Gx)M: 

_L_^(wt"w'"), 



is the characteristic function of the random variable S, Cs{s), 
sampled in s = — j27rfc. In particular when (5 is a zero 
mean Gaussian random variable with variance ag, we have: 
C = cxp(cr|WV2) and (C)fcfc = ciip{-2TT^Pal), k = 
-M, . . . , Af. 

b) Computation of k[GI.Gx]-' 
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which does not depend on k. Thus, diag(Z) 
_i_Tr { Wt " W" } I, and 
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where Z = GtWt"W™Gx. Now, 
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Finally, 



E ^V^diag(Z) 
n!m! 
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] = Ai("V^'"UZ-diag(Z))+/^("+™)diag(Z) 
^(»)^(™)Z + (^(«+™) - diag(Z) 
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*Let A = diag(ai, . . . , a,i) be a diagonal n X n matrix. The exponential 
of A, denoted by exp(A), is the diagonal matrix whose elements are 
[exp(ai), . . . ,Gxp(a„)]. 
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where 
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Therefore, 



Concluding, 
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E [Gt G.] = Gt ctCG. + (1 - I 
is real, 

[GtGJ = GtC^Gi +(l- Ili^] I 



and, if C is real. 



